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ABSTRACT 

We study low-amplitude crustal oscillations of slowly rotating relativistic stars con- 
sisting of a central fluid core and an outer thin solid crust. We estimate the effect of 
rotation on the torsional toroidal modes and on the intcrfacial and shear spheroidal 
modes. The results compared against the Newtonian ones for wide range of neutron 
star models and equations of state. 
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INTRODUCTION 



Crustal oscillations of nonrotating nonmagnetic neutron stars have been studied in New t onian theory ( Hjmsen fc Cioffil 198pl : 



McDermott et al. 1985 



General Relativity (jSchumaker fc Thorndll98. 



McDermott. Van Horn fc HanserJll988l: [&t7ohmaver et al.lll99ll: iBastrukov et al.1 20071) as well as i n 
i lFinJl990l ; lLeins!ll994l : lYoshida fc Leelbood: ISamuelsson fc Andersson1l2007l) 



These studies have been gradually ext ended in order to include the effects of rotation (jstrohmavei 1991 : iLee fc Strohmaver 

1996| ; lLeel20073;lyavoulidis" et alj 20071) or strong magnetic fields JCarroll et al Jl986l:lDuncanll998l:lMessios. Papadopoulos fc Stergioulal] 



2001 



Pirol 



2005; Lee 2007a 



7 al: ISotani. Kokkotas fc Stergioulasll2007al : lLee!l2007bl : 



Sotani. Colaiuda fc Kokkotas! 120071 ). Most of 



them have been focused on torsional toroidal modes of oscillation (designated as it n ) and only a few studies have dealt with 
interfacial and shear (is n ) spheroidal modes. There are also some recent studies of the perturbations of pu rely elastic 
stars in the general relativistic framework (jKarlovini. Samuelsson fc Zarroual2004l ; iKarlovini fc Samuelssonll2007l ). 

After the di scovery of high-frequency quasi-periodic oscillations (QPOs) in the tails of giant flares from soft gamma-ray re 



peate rs (SGRs) (jf3arat et al.lll983l : llsrael et al.ll2005l : lstrohmaver fc Wattsl2005l : lwatts fc Strohmaverl2006l ; lstrohmaver fc Watt's! 
2006). special attention has been drawn to torsional toroidal oscillations of neutron star crusts. Most of the QPOs have been 
observed at frequencies between 18 and 155 H z, although there have been a few at higher frequencies, for example at 625, 
1840 and possibly at 718 Hz for SGR 1806- 20 (jWatts fc Strohmaverll2007l ). Fundamental torsional modes could account fo r 
many of the low frequencies ( Duncan 19981 ) while radial overtones could account for some of the higher ones ( Piro 2005 ). 
However, the torsional-mode interpretation had two drawbacks: it could not explain the very low observed frequency at 18 
Hz and it could not explain observed pairs of frequencies as those at 26 and 30 Hz or those at 625 and 718 Hz. 

It has been soon realized that the coupling between the crust and the magnetic field would play a key-role in our 
attempt to explain the observed QPOs. In fact, this coupling would favour the existence of global magnetoelastic modes of 
oscillation rather than pure elastic modes confined in the crust I Levin 2006h . Moreover, it has been argued that these modes 
should decay on a short timescale because of the presence of a magnetohydrodynamical continuum in the c o re and only 
specific QPOs could be long-lived dGlampedakis. Samuelsson fc Andersson 20061 ; Sotani et al. 20061 : Lee 2007a ; Levin 2007 ; 
Sotani. Kokkotas fc Stergioulasll2007bl ). 

After a catastrophic reconfiguration of the stellar magnetic field, axial-type torsional and Alfven oscillations should be 
the most easily excited as polar-type oscillations would have to overcome strong restoring forces. However, the presence of 
a magnetic field would inevitably couple these oscillations with polar-type ones characterized by interfacial (i-), shear (s- 
), pressure-restored (f- and p-), gravity (g-) and polar-type magnetoelastic modes; stellar rotation should do more or less 
the same. These polar-type modes would involve density variations and could be relevant for gravitational-wave emission 
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( Abbott et al.l 2007, LI GO scientific collaboration). The po ssible association of polar- type modes with the observed QPOs 
has been hinted by Pirol (|2005h ; IWatts fc Strohmaverl (|2006h but further investigation is definitely needed. 

In Section[2]we derive the equations that describe linear spheroidal and toroidal oscillations of slowly rotating relativistic 
stars possessing a crust. In most of them, a parameter e indicates the presence of a relati vistic term. Therefore e — > becomes 
our s tandard tool to recover the already known Newtonian equations for these oscillations (jStrohmave rll99ll :l iLee fc Strohmaverl 
19961). Furthermore, the nonrotating parts of these relativistic equations are easily comparable with the equations given in 



Yoshida < 



Leel J2002S). In Section H we derive the perturbation equations for the solid-crust region (|2.1[) . we review those for 
the fluid-core region (|2.2[) and we present the necessary physical boundary and jump conditions (|2,3fl . In Section[3]we present 



our numerical results for a set of neutron star models with different equations of state (EoS) and different bulk properties 
(13.11 1. first revisiting torsional toroidal modes (|3.2p and then focusing on interfacial and shear spheroidal modes (|3.3|l . In 
Section U we summarize and discuss our results. 



2 FORMULATION 



We consider a slowly rotating relativistic, strain free, star described by the metric: 



2i- / ji2 . 2A t 2 . 2 t/i2 . 2-2 

-e at + e dr +r do +r sin 



2uir sin 8dtd<j>, 



(1) 



where v, X and uj are functions of the radial coordinate r. These functions are solutions of the Tolman-Oppenheimer-Volkoff 
(TOV) equations: 



dv 
dr 
dX 
dr 



_2M (r 



4nrep + 



M(r) 



dM (r) 2 

- ± - L = VKT p, 

dr 

dp dv 
_ = _ (p + ep )_ 



ee 2X | 4nrp - 



M(r) 



and of one more equation describing the dragging of the inertial frames of reference ( Hartlei 1967 ) 



d uj ( dv dX 
dr 2 \ dr dr 



duj 
dr 



16-7re (p + p) uj — 0, 



(2) 
(3) 
(4) 

(5) 



where w := O — uj, p and p are the energy density and the pressure, respectively, M (r) is the mass inside radius r and Q is 
the stellar rotational frequency. In our general-relativistic approach, the parameter e equals to 1. Its presence aims to point 
out the Newtonian limit where e — > 0. Then, according to equations ©-(HI), X — > 0, dv/dr — > M (r) /r 2 , dp/dr — > —pdv/dr 
and dX/dr — > 0. Furthermore, in the Newtonian limit, uj — > or, equivalently, zu —* Q. and dzo/dr — > 0. 

The pulsation equations come from the linear perturbation of the energy-momentum conservation law 8 (V^T Q ^) = 
where: 



(p + p) u a U{3 + pg afj - 2pLS a 



(6) 



The components of the perturbed four-velocity are given by the relation 8u a = £ u £ a where £ a is the displacement vector 
while the shear tensor S a /3 is given by the relation a a g = £ u S a g where a a g is the rate of shear tensor. Namely a a p is the Lie 
derivative of the shear tensor along the world lines ( Carter fc Quintana 1972} ) and is calculated by the equation: 



where P t 



= 2 



is the projection tensor: 



Pa/3 = g a f3 + U a Ug. 



(7) 



(8) 



It is implied that we work in the Cowling approximation as we neglect the perturbed Einstein equations SG a /3 = 9>n8T a g and, 
additionally, set all metric perturbations equal to zero in our equations. The Cowling approximation is typically very good 
for toroidal type of oscillations i.e. t-modes, r-modes etc, while typically for polar type of perturbations the error can be of 
the order of 10-20% especially for the fundamental pressure mode the /-mode. For, the other type of modes, like p-modes, 
higher p- modes and s-modes one does not expect deviations larger than 3-5%. 

We choose to work in a corotating reference frame where 8u a — e~ v d^ a / 'dt, a a g — e~ v dS a fS / dt. Using equations (|6])-(|8|), 
the energy-momentum conservation law yields three second-order partial differential equations for the three components of 
the displacement vector £ r = £ r (t, r, 8, (j> \ £® = £ g (t, r,8,(f>) and £^ = ^ (t, r, 8, 4>). These equations are the relativistic version 
of the Newtonian equations (35)-(39) of IStrohmaverl ( 19911 ) and their explicit form is given in the Appendix lAl 



When we focus on spheroidal modes, we can write for the displacement vector: 
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e = 



3 1 d 
oO sm 2 6 o<p 



while when we focus on toroidal modes, we can write: 



e 



0,T 



J. d_ 

sin 8 d<j) ' 



-T- 







(9) 



(10) 



Working in the slow-rotation approximation (jKoiimal 1 19921 ; IStavridis fc Kokkota s 2005), we eventually obtain three second- 
order ordinary differential equations for the functions S — S (r) , H — H (r) and T = T (r). The equations for S and H are 
quite le ngthy and we do not sho w them here. However, we mention that they are the relativistic version of the equations 
used in Lee fc Strohmavei ( 19961 ). On the other hand, the equation for the toroidal radial function T, being shorter, reads 
( Vavoulidis et alJl2007T h 



2rn 2 2e(v-\) 
a 1 = v B e 



St 

dr 2 



4 dv dX 1 dfi\ dT 
r dr dr /x dr J dr 







— Imaw 





T. 



(11) 



where v'j, := /i/ (p + ep) is the speed of shear waves and A := L (£ + 1). 

We then expand the frequency and the displacement functions in power series: 



o- = £ rfcTj = a + Wl + O [rf) , ? = £ rfC' j = t° + vC' 1 + O (r? 2 ) 

3=0 3=0 



(12) 



where rj is an auxiliary expansion parameter and we insert these expansions into the main perturbation equations, e.g. equation 
(|11[) above. In this way we split each equation in two, one zeroth-order in SI which is taken collecting all 77° terms and one 
first-order in 17 which is taken collecting all rj 1 terms. For example, equation (111[) , to zeroth order in rotation, gives the 
equation: 



2rrn0 

a T 



2 2e(v-\) 



d 2 T" 



dr 2 



+ 



4 dv 
r dr 



d\ + 1 dfi\ dT 
dr /j, dr 



dr 



,A-2 t q 



(13) 



which, when supplied with the appropriate boundary conditions, determines the zeroth-order eigenfrequency, cro, and the 
zeroth-order eigenf unction, T°. Then, to first order in rotation, equation (|ll|l gives another equation: 



2 r7~i 1 o 

<7q1 — L(7qO\± 



d T 
dr 2 



dX t 1 dp\ dT 1 
dr 



dr dr ~^ fi dr 



A-2 



T 1 



1 - 



A 



(14) 



which determines the first-order rotational corrections of the eigenfrequency, a\ , and of the eigenfunction, T 1 , using the known 
zeroth-order quantities cro and T , already calculated by equation (|13p . 



2.1 Solid crust 

To describe the oscillations in the solid-crust region, we use the following functions: 

4 = s>, 



r 2 
r--ai 



-2eA dH ] j 

e r — ■ h b 

dr 



-2,A dP 

aie r . 

dr 



(15) 
(16) 
(17) 
(18) 
(19) 
(20) 



This set of functions is the set (31)- (36) of lYoshida fc Leel (|2002h and it reduces to the set of functions (53), (16), (15) of 
Lee fc Strohmaverl (|l996l ) in the Newtonian limit, 



0. 



Using these functions, we can recast our three main perturbation equations into a system of six first-order ordinary 
differential equations where four of them are describing the spheroidal perturbations while the other two are describing the 
toroidal perturbations. To zeroth order in Q (j = 0), this system has the form: 

r d A = -(i + 2^+eU 2 ) Z 1 + ±-z 2 + ^e(£+l)zl (21) 



03 



a 3 
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,d4 
dr 



dr 



+ 



-3 - eU 2 + Ui - e 2£A cio- 2 ) Vi + 4^- (3a 2 + 2ai)| z? + (V 2 - 4^ij z% 
Vi - 2ai + 2g j | * (/+ 1) 4 + e 2sA 7 (* + 1) z\, 

„2eA 
2eA , e 
— e 2! H Z 4 , 

ai 



$i = - (-Vi + 6r^i") 2 ? - ^2 2 ° - ( Cl a 2 Vi + 2ai - 2^1 (a 2 + a 3 )i (< + 1)1 4 ~ (3 + etfc - V 2 ) 24°, 



dr 
cir 
cir 



2eA 
z 6> 



(22) 
(23) 
(24) 
(25) 
(26) 



- {cia 2 Vi - 01 (£ - 1) (£+ 2)} z\ - (3 + et/ 2 - V 2 ) z° 6 . 

This set of equations is the set (25V(30) oflYoshida fc Lee! (|2002h and, in the Newtonian limit, reduces to the sets of equations 
(54)-(57) and (68)-(69) of ILee fc Strohmaverl Jl996h ~ The first four equations give the zeroth-order eigenfrequency ao and 
eigenfunctions Z\, z 2 , 23, z\ for the interfacial and shear spheroidal oscillations while the last two equations give ao, z% and z% 
for the to rsional toroidal o scillat ions. 

As in lYoshida fc Leel |2002h . the functions that appear in equations (|21|l 



cti = — , a.2 = T— -qi, «3 = r+-ai, 
p 3 3 

p + ep df t t- P dv 

Vi = ^ r — , V 2 = -r — , 

p dr p dr 



Ui 



1 d / di^ 



dr V dr 



U 2 =er 



dX 
dr ' 



du\ 1 



Af _ 2e „ . 
= Jp re [dr-) 

In the Newtonian limit, e — > and dv/dr - 
Therefore, the functions (1281)- (1301) reduce to: 



(27) 
(28) 
(29) 
(30) 

where g = M (r) /r 2 is the gravitational acceleration with dp/dr — ► — pg. 



p dlnp 
Vi, Vi> -> -rp = -— — 
p dmr 



= V, 



d (M(r)\ d\aM(r) 



dlnr 



ci 



M 



rg 



-'=( 



1 = 17-1, 



r \ 3 JVif 



i?7 M(r) 



(31) 
(32) 
(33) 



The fu nctions V, U and Ci are well- known from the Newtonian theory, see e. g. the set of relations (35) of ILee fc Strohmaver 
I 1996h or the set of relations (15) of iMcDermott. Van Horn fc Hansen! (jl988T ). 
To first order in SI (j = 1), our system of equations has the form: 



dr 
, dz 2 
dr 



1 + 2^ + ef/a^ zl + —z\ + + 1) 4, 



(34) 



-3 - eU 2 + 7/i - e 2£A Cl( j 2 Vi + 4— (3a 2 + 2ai) ^ 21 + V 2 - 4^- 2^ 
Vi - 2ai (l + 2^ j £ (7 + 1) 4 + e 2eA ^ (7 + 1) z] 
+ j— 2e 2eA ci<T a"i Vi + e»4j 2? + {2mci<7oroVi + eS} 23 + eC^, 



,d4 
dr 

_dz:[ 
dr 



, dgs 
dr 

dzg 



2 £ A 

2e\ 1 . e 1 

-e 2! H z 4 , 

ai 



+ 



-Vi + 6r— *} - — z\ - \ Cl a 2 Vi + 2 ai -2—{a 2 +a 3 )l(£ + l)\ z\ - (3 + et/ 2 - V 2 ) z\ 

«3 / Q3 ( Q!3 J 

f 2mciCT roVi \ , c , f T/ , 2mciao^Vi 1 

j + eD | Zl + £ ^ 2 + |- 2c ^o^^ + H£+1) + 2 3, 



(35) 
(36) 

(37) 
(38) 



= - { Cl a 2 Vi - ai (£- 1) (i + 2)} 4 -(3 + eU 2 -V 2 )z\ 
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where: 

A = 

B = 

C = 

V = 

£ = 

T = 

Q = 



2cio- Vi 



-o-i + 



1(1- 



+ eg\zl 



(qi 
P 
,P 

-2v+2\ I 013 



CX2) r maovj, 

\ 3 dvj 
+ 1 + Q2 I r ma — 



— 1 ) r maozu, 

ai 



£ + l-ai 
P 



r mao dvo 
£(£+1)17 



i (r 



CK2 



dv 2r dp , . 
r— + --- r- +2(ai - a 2 ) 



dr 3 p dr 



r moazo 



dv r dp 

1 — ai r— H — + 2 Qi — a 2 + «3 + 

ar p dr 



OL\OL2 
«3 



r mffon 



^+1) 



-2k / «i 



_L _ i 

a 3 



1(1- 



— 4ai + 2a 3 — an + 



a :1 



« (-£ + !) 



r 771(70 

^ + 1) ' 



'2a! -!)(£ + 2) 



r maozu 



(39) 

(40) 
(41) 

(42) 

(43) 

(44) 

(45) 

(46) 



£(l+l)' 

From equations (1341) - (|39[ ) we get the first-order rotational correction of the eigenfrequency, <7i , and the first-order rotational 
corrections of the eigenfunctions, z\ (i = 1 . . . 6). Again, the first four of them refer to the interfacial and shear spheroidal 
modes while the last two refer to the torsional toroidal modes. Bars over oyj, o\ and vj indicate dimensionless quantities, scaled 

(To I \J M/R 3 . This set of equations reduces to the sets of equations (58)-(61) and (70)-(71) of 



by \/M/R?\ for example oyj 



Lee fc Strohmaverl (jl996l ) in the Newtonian limit. 



2.2 Fluid core 

To describe the oscillations in the fluid-core region, we use the following functions: 

y\ = S\ 



dv 
dr 



dp 3 



p + ep 



+o 5p J 
grp' 



The zeroth and first-order spheroidal radial functions H° and H 1 are given by the relations 
H° = 

H 1 = 



-9 ' 



f 2mzu 



V 



(TO 



y? 



^ + 1) 



+ 



2mvj 



mixj 



+ e 



£(£ + 1)' 



— 9 ' 



(47) 
(48) 

(49) 
(50) 



1(1 + 1) 

cf. rel ations (44) of lYoshida fc Lee (|2002l ) for nonrotating relativistic stars and relations (81),(86)-(87) of iLee fc Strohmaverl 
1 19961 ) 'or rotating Newtonian stars. 

To zeroth order in Q (j = 0), the system of equations ()Alf) - (|A3() with p = can be recasted into a system of two first-order 
ordinary differential equations for the functions y( and y®- 

r M 

dr 
dy\ 



Vi 



2/2 ! 



(e 2£A cia 2 + rA r ) yl - (C/i + rA r ) y° 2 , 

where A r is the relativistic Schwarzschild discriminant: 
1 dp 1 dp 



(51) 
(52) 



.4 



p + ep dr Yp dr 

Similarly, to first order in fi (j = 1), our system for the functions y\ and y\ is the following: 



(53) 



,dyl 
dr 



= (e^c^l + rA^yl-tU.+rA^yl 



3-^+eU 2 



+ 



2mzu 



■eH\y 



U 

r 

2mm 



f 1 - I y 



' + 1) 



1 



(Jo 



£(£ + l)^ 



+ 2el 



JtL 

— 9 ' 

ClOo 



(54) 
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n 2eA - - 

+ 2e aaoaiyi - 



2mzo 



eH}y u 2 



where: 
H = 
1 = 

In the Newtonian limit s 



dzu 
dr 



dr 



rm 



-2l/ 2 



(55) 

(56) 
(57) 



(84)- (85) of 
(42)- (43) of 



Vj — > V, — ► Q and equations (f5"l"j) - (f52)) and ([54l) - (f55)l reduce to equations (82)-(83) and 
Lee fc Strohmaverl dl996h . respectively. In the relativistic nonrotating limit, — > 0, we are reduced to equations 



Yoshida fe Led l|2002h 



2.3 Boundary and normalization conditions 

At the stellar center (r = 0), the eigenfunctions j/^ and j/j must be regular. By expanding them in appropriate series, t/i 
2~2o a i,nT n and carrying out some algebraic manipulations, we find: 

cialv{ - £y 3 2 + F 3 



0. 



where: 

F° = 0, 
i 2muj I a\ 1 



-2 o 

cicr y 1 . 



(58) 

(59) 
(60) 



At the stellar surface (r = R), the Lagrangian perturbation of the pressure must vanish (Ap = 0). This eventually leads 
to a simple relation between y\ and y\, i.e.: 



o. 



At the fluid-solid interfaces, we require continuity of the tractions. This means the jump conditions: 

7 i 
z{ = VI, 

4 = Vi(y(-yi), 

4 = o, 

for the spheroidal modes, and: 



for the toroidal modes. 

Finally, we normalize our zeroth and first-order eigenfunctions by imposing the conditions: 

y°i = l, 
yl = o, 

respectively, at the stellar surface. 



(61) 

(62) 
(63) 
(64) 

(65) 



(66) 
(67) 



3 RESULTS 



3.1 Neutron star models 



We w ork with a set of 34 realistic neutron star models: we combi ne EoS A ( Paiidliari£andell97l|) , WFF3 (jWiringa. Fiks fc Fabrocini 
1988h. APR dAkmal. Pandharipande fc Ravenhalllll998ri or L (jPandharipande fc Smithlll975h for the fluid core with EoS DH 
( Douchin fc Haensell 200ll ) or NV jNegele fc Vautherin 19731 ) for the solid crust and, for each combination, we construct a 
sequence of models, beginning from a low-mass model of 1.4Mq and reaching, in increments of O.2M0, the maximum-mass 
model allowed by that EoS. That maximum-mass model is > 1.6, 1.8, 2.2 an d 2.6M W for EoS A, WFF3, APR and L, respec- 
tively. Further details about this set of neutron star models can be found in Sotani. Kokkotas fc Stergioulasl (2007a) and in 
references therein. Here, we add the following pair of fitting formulas: 



In (Ar/R) ~ -7.95 (±0.04) M/R - 1.28 (±0.01) , for the DH crustal EoS, 
In (Ar/R) ~ -7.84 (±0.03) M/R - 1.04 (±0.01) , for the NV crustal EoS, 



(68) 
(69) 
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which relate stellar compactness M/R with relative crust thickness Ar/R. It is obvious that this formulae are valid for typical 
neutron stars with a fluid core and a crust. In the limit where M/R — > the star "solidifies" or else the crust equation of state 
is valid for the whole star and the above formulae do not reach the correct l i mit w hich is just 1. We should point out that a 
semi-analytic formula has already been derived by Samuelsson fc Anderssonl 1 20071 ) . this formula is valid for a wide range of 
stellar compactness and has t he correct M/R — > limit. In thi s formula [eq (B6)] there is a free parameter a which should be 
fixed by the equation of state. Samuelsson fc Anderssonl ( 2007 ). for the crust EoS DH, find a ~ 0.019 for F = 4/3 polytropes 
and a ~ 0.023 for the set of realistic equations of state for the fluid core that they have used. For our stellar models we find 
a « 0.020 for the crust EoS DH and a = 0.026 for the crust EoS NV. 



3.2 Torsional modes 

The Newtonian fir st-order rotational corrections of the eigenfrequencies of the torsional modes are shown in the formula of 



Strohmaverl l|l99lt l 



a = a + o-i = <T + mS!C liNew = a + Jffpfi > ( 70 ) 

which indicates that C-y N ew '■— ai/mSl takes the value 1/ [1(1+ 1)] independently of the structure of the star. For example, 
Cy N ew = 0.167 for £ — 2, ]\f ew = 0.083 for £ = 3, N ew = 0.05 for £ = 4 and so on. Moreover, the Newtonian first-order 



rota tional corrections of the eigenfunctions of the torsional modes can be set equal to zero, T 1 — 0, cf. IStrohmaverl (|199lf ) 



and Lee fc Strohmaver ( 19961 ) . The relativistic results can be drawn from Tables Q] and [5] and from Figure [T] The relativistic 



first-order rotational corrections of the eigenfrequencies of the torsional modes, as measured from a rotating observer, are 
shown in the formula: 

a — ao + ai — ao + mfiCS re i . (71) 

Then, the rotating observer compares the Newtonian and the relativistic first-order rotational corrections defining rcr := 
1 — Cy le [/Ci N ew - On the other hand, an inertial observer measures: 

a = ao - mO, + mnC 1)New = ao + mfiC^ New , (72) 
a = (To — mil ± mflCy re [ = ao ± mQC'^ re i, (73) 

in the Newtonian and in the relativistic case, respectively, and he defines rq := 1 — C'^ m \/C'i New' ^ * s eas y t° check that 
rcr and rcj obey the simple relation: 

(74) 



1 / <7 l,New _ 1 

In Table [T] we list the eigenfrequencies ao and the relativistic first-order rotational corrections C-^ r „j for the fundamental 
I = 2, 3 and 4 torsional modes (2*0, 3^0, 4^o) for our set of 34 stellar models. The Newtonian first-order rotational corrections 
N ew are equal to 1/ [£(£+ 1)] for all stellar models according to Strohmayer's formula (|70p . As expected, for each EoS, 
higher compactnesses M/R result in higher relativistic corrections rcr. For the chosen set of stellar models, where 0.14 < 
M/R < 0.28, the relativistic corrections rcr vary from 10% to 30% approximately. According to equation (|74|) and for i = 2 
{pl New = 0-1(57^ , an inertial observer measures relativistic corrections rcj = — rcr/5, that is from -2% to -6% approximately. 



These results are in good agreement with earlier investigations (|Vavoulidis et al.l 120071 ). In Table [2] we turn our attention to 



higher (n = 1, 2, 3) torsional overtones. Dipole (£ = 1) modes are now allowed by the angular-momentum cons ervation law and 



we additionally know that these higher overt ones are quite insensitive to the azimuthal harmonic index £ IjHansen fc Cioffi 



1980l : iMcDermott, Van Horn fc Hansenlll988l ) . Table [2] refers to £ = 1 but results for I > 2 are very much alike. In Appendix 



[Bl we briefly show that the relativistic first-order rotational corrections of the eigenfrequencies of the torsional modes can be 
estimated also by the integral formula: 

f,f w [l/A + ev 2 s (1 - 2/A)l (T ) 2 dr 
a 1 = mtlC l rel = 1 ' ° K ' ^ ' • (75) 



The numerical evaluation of formula (|75[) and the numerical solution of the previously described eigenvalue problem yield 
essentially the same results for a\. 

Finally, in Figure [T] we show the relativistic corrections rcr versus compactnesses M/R for our set of 34 stellar models, 
based on the results of Table [T] Clearly, models with the NV crustal EoS are shifted towards higher rcr in respect with models 
with the DH crustal EoS; linear fits reveal that: 

rc r (%) ~ 108.9 (±3.0) M/R - 3.5 (±0.7) , for the DH crustal EoS, (76) 
rc r (%) =± 108.9 (±4.9) M/R - 2.5 (±1.1) , for the NV crustal EoS, (77) 
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Table 1. Eigonfrequencies (otj), Newtonian first-order rotational corrections y^i N CW J an< ^ rclativistic first-order rotational corrections 
(C\ re l) of the £ = 2,3 and 4, n = torsional modes (2<o> 3*0! 4*o) f° r various stellar models. As expected, for each EoS, higher com- 
pactnesses (M/R) result in higher rclativistic corrections ^rcr := 1 — ^ e \/C\ ]\J CW )- The Newtonian first-order rotational corrections 
are equal to !/[£(£ + 1)] for all stellar models. 



Model 


M/R 
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c 
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vo (Hz) 
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and these can be combined with the fitting formulas (|68p ~ (|69p of the previous Subsection. 



3.3 Interfacial and shear modes 



The Newtonian first-order rotational corrections of the eigenfrequencies of the interfacial and shear modes can not be given by 
an analytic formula like Strohmaver's formula (|70[1 . However, these corrections are given by the integral formula (|Unno et al 



19891 ; IStrohmaverlll99ll ; iLee fc Strohmaverll 19961 ) 



C 



l,New 



So* Pff + 



(c h ) 2 ] 



pr 2 dr 



(78) 



where C = rS°, $ h = rH° are the zeroth-order radial eigenfunctions of the displacement vector (compare with equation [9}. 
We clearly see that the Newtonian first-order rotational corrections of interfacial and shear spheroidal modes do depend on 
the background stellar model whereas torsional toroidal modes were proved to be independent. 

In Table[3]we list values of Ci for the I — 1,2 and 3 interfacial modes for our set of 34 stellar models. ]\f ew is calculated 
by evaluating the integral formula (|78|) while Ci re [ is calculated by solving the eigenvalue problem (|34j) - f|3T[> . (|54[1 - (|55|) along 
with the appropriate boundary and jump conditions. In both cases, we first need to solve the zeroth-order eigenvalue broblem 
in order to determine the zeroth-order eigenfrequency and the zeroth-order eigenfunctions. 
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Table 2. Eigenfrequencies (cro), Newtonian first-order rotational corrections \ C\ N CW J ana - relativistic first-order rotational corrections 
(C\ re l) of the i = 1, n = 1, 2 and 3 torsional modes (iti, 1*2, 1*3) for various stellar models. As expected, for each EoS, higher com- 
pactnesses (M/R) result in higher relativistic corrections ^rcr := 1 — New) - '^ ne Newtonian first-order rotational corrections 
are equal to 1 /2 for all stellar models as t = 1 . 



Model 


M/R 


era (Hz) 


n = 1 
^-T,New 


w l,rel 


an (Hz) 


n = 2 
^-T,New 


w l,rel 


cro (Hz) 


n = 3 
^-T,New 


C1 1 
^l,rel 


rcr (%) 


A+DH14 
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0.5 
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0.5 
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0.5 
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18.89 


A+DHi 6 
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0.5 
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0.5 


0.379 
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0.5 
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24.30 
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0.191 
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0.5 


0.419 
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0.5 


0.420 
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0.5 
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16.11 


WFF3+DHi 6 


0.223 


1099.0 


0.5 


0.402 


1834.0 


0.5 


0.403 


2531.8 


0.5 


0.403 


19.50 


WFF3+DHis 


0.265 


1363.8 


0.5 


0.378 


2275.3 


0.5 


0.378 


3144.0 


0.5 


0.378 


24.48 


APR+DH14 


0.171 


760.0 


0.5 


0.429 


1268.9 


0.5 


0.429 


1748.2 


0.5 


0.430 


14.12 


APR+DHie 


0.195 
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0.5 


0.417 
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0.5 


0.417 


1976.6 


0.5 


0.417 


16.69 


APR+DHi 8 


0.221 
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0.5 


0.402 
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0.5 


0.402 
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0.5 


0.403 


19.59 
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1080.9 


0.5 


0.385 


1803.4 


0.5 


0.385 


2490.8 


0.5 


0.386 


22.95 


APR+DH22 


0.279 


1234.7 


0.5 


0.364 


2060.1 


0.5 


0.364 


2846.4 


0.5 


0.365 


27.17 


L+DH14 


0.141 


529.2 


0.5 


0.441 


884.2 


0.5 


0.441 


\2\7A 


0.5 


0.442 


11.83 


L+DHi 6 


0.160 


585.3 


0.5 


0.432 


977.7 


0.5 


0.432 


1347.9 


0.5 


0.433 


13.61 
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0.179 
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0.5 
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0.5 
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0.5 
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0.5 
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0.5 
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0.5 
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0.5 
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0.5 
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0.5 
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0.5 
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0.5 
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1365.3 
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0.5 
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L+NV14 


0.152 
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0.5 


0.428 


858.5 


0.5 


0.429 


1212.3 


0.5 


0.430 


14.35 


L+NVie 


0.171 


524.0 


0.5 


0.419 


931.1 


0.5 


0.419 


1316.6 


0.5 


0.420 


16.21 


L+NVig 


0.190 
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0.5 


0.409 
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0.5 


0.409 


1426.1 


0.5 


0.410 


18.16 


L+NV20 


0.210 


614.5 


0.5 


0.399 
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0.5 


0.399 


1545.8 


0.5 


0.400 


20.25 


L+NV22 


0.230 


666.7 


0.5 


0.387 


1183.9 


0.5 


0.387 
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0.5 
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22.55 


L+NV24 


0.253 
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0.5 
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1294.0 


0.5 


0.374 


1835.2 


0.5 
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25.23 


L+NV 26 
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823.8 


0.5 


0.356 


1462.1 


0.5 


0.356 


2074.1 


0.5 


0.357 


28.74 



In Table |4] we focus on the shear modes. These eigenmodes are quite insensitive to the value of I, thus we fix I = 1 and 
list the Newtonian and the relativistic first-order rotational corrections of the first three members of this family of modes 
is 2 , 1S3). We note that the eigenfrequencies of the shear modes are almost equal to those of the higher overtones of the 
torsional modes. This is expected as s-modes and t-modes represent different polarizations of transverse elastic waves. 



4 DISCUSSION 

In this paper, we have calculated the first-order rotational corrections of the torsional, interfacial and shear modes of realistic 
neutron stars with crusts within the framework of General Relativity. Our results agree, qualitatively, with previous Newtonian 
results but differ, quantitatively, a few per cent. An interesting observation is that, although the low-frequency torsional modes 
could not explain all of the 18, 26 and 30 Hz frequencies observed in SGRs, there is a good chance that one or two of them 
could be explained by low-frequency interfacial modes as it can be seen from Table [3] This would be a novel scenario 
for gravitational-wave detection since it would imply the excitation of polar-type modes via their coupling with torsional 
oscillations. 

Rotation could drive many oscillation modes unstable through the Chandrasekhar-Friedman-Schutz (CFS) mechanism. 
For an observer rotating with the star, the instability would set in when the eigenfrequency of a mode oy would be equal to mil; 
at the same time, the eigenfrequency of the same mode for a distant inertial observer en would be zero. For the torsional and the 
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Figure 1. Relativistic correction as measured from a rotating observer (rcr) versus compactness [M/R) for the fundamental torsional 
modes ((to) for various EoS. 



inter facial crustal modes this would happen at moderate stellar rotational frequencies (jYoshida fc Ledl200ll ; IVavoulidis et al 
2007|). For example, in Figure [5] which refers to the L+DH14 EoS, we show that this would happen at Q ~ 12.5Hz for the 
torsional 2^0 mode and at £} ~ 19Hz for the interfacial i% mode. We have adopted the slo w-rotation approximation b ut this 



proves to be sufficient in most cases. For example, SGRs rotate with periods of ~5-8 sees ( Woods fc Thompson 20061 ) which 
are extremely high compared to their Kepler periods (>lms). On the other hand, for such slow rotators, rotational effects are 
quite small; it would be a big surprise if they could be measured from current x-ray detectors. Furthermore, CFS instabilities 
of crustal modes could not be developed in those objects as they would require stellar rotational frequencies of several tenths 
to hundreds of Hz to set in. However, for newly-born neutron stars things may be quite different. Newly-born neutron stars 
are expected to rotate much faster, with frequencies from tenths to hundreds of Hz, and therefore they could become CFS 
unstable soon after they would form their crusts. In such a case, the slow-rotation approximation would be still valid while 
the effects of rotation in the spectrum would become pronounced. 
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Table 3. Eigonfrequencies (cro), Newtonian first-order rotational corrections (C^ ]\J CW J and relativistic first-order rotational corrections 
(C\ re l) of the t = 1, 2 and 3 interfacial modes 2h 3*) for various stellar models. 
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Table 4. Eigenfrequencies (<xo), Newtonian first-order rotational corrections (c^ ]\J CW ^ an< ^ rclativistic first-order rotational corrections 
(C\ re l) of the t = 1, n = 1,2 and 3 shear modes (isi, 1S2, 1S3) for various stellar models. 
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APPENDIX A: RELATIVISTIC VERSION OF STROHMAYER'S (1991) EQUATIONS 



Here we give the explicit form of the three second-order partial differential equations for the three components of the displace- 
ment vector £ r = £ r (t, r, 9, <j>), (f = £ e (t, r, 9, <j>) and = £f (t, r, 9, <j>). In order to be able to compare easily with the Newto- 
nian equations (35)- (39) of Strohmavei ( 199ll ). we define the displacement vector as £ l 
and we get: 
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APPENDIX B: PROOF OF FORMULA ([75)1 



Here we apply a mathematical technique that allows us to estimate the relativistic first-order rotational corrections of the 
eigenfrequencies of the torsi onal modes through a simp le integral formula. The same technique has been applied in fluid slowly 
rotating relativistic stars bv lYoshida fc Koiimal (|l997l ). 

Multiplying equation (|14|l with the conjugate function T°* and integrating over the whole star we get: 
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which implies that the first term in equation (|B1[) annuls, leaving the following integral formula for the corrections o\\ 
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In the Newtonian limit (e — > 0, w — > fl), this formula becomes: 

Q ^ A) = WTTy (B4) 

as expected. 



